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We present a measurement of the IV-boson-pair production cross section in pp collisions at 1.96 
TeV center-of-mass energy and the first measurement of the differential cross section as a function 
of jet multiplicity and leading-jet energy. The W + W~ cross section is measured in the final state 
comprising two charged leptons and neutrinos, where either charged lepton can be an electron or 
a muon. Using data collected by the CDF experiment corresponding to 9.7 fb -1 of integrated 
luminosity, a total of 3027 collision events consistent with W + W~ production are observed with 
an estimated background contribution of 1790 ± 190 events. The measured total cross section is 
o(pp —> W + W~) = 14.0 ± 0.6 (stat)lj Q (syst) ± 0.8 (lumi) pb, consistent with the standard model 
prediction. 
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The measurement of W-boson-pair production 
is an important test of the electroweak gauge 
sector of the standard model (SM) of particle 
physics GLi- Extending the measurement to 
include the properties of jets produced in associ¬ 
ation with the boson pair is an interesting test of 
quantum chromodynamics (QCD) that has not 
been performed before in events with multiple 
heavy gauge bosons. This process is a signifi¬ 
cant background for Higgs-boson studies in the 
W + W~ (WW) decay mode at particle colliders, 
where the use of jet vetoes and jet counting is 
an essential element of the analysis technique [3|— 
Jfji]. Furthermore, a final state consisting of two 
massive gauge bosons and two jets is typical of 
vector-boson scattering, a process sensitive to 
many non-SM contributions to the electroweak- 
symmetry-breaking sector Q. With the first ev¬ 
idence of vector-boson scattering at the Large 
Hadron Collider (LHC) @, i , it is essential to 
verify and tune the simulation tools used to de¬ 
scribe electroweak processes produced in associ¬ 
ation with jets. 

This communication reports measurements of 
the W + W~ production and differential cross sec¬ 
tions as a function of jet multiplicity and jet en¬ 
ergy in a final state consisting of events with two 
oppositely charged leptons, where each lepton is 
identified as either an electron or a muon, and 
an imbalance in the total event transverse mo¬ 
mentum (transverse energy), due to the presence 
of neutrinos. The differential measurements of 
jet multiplicity and jet energy are unfolded to 
hadronic-jet level to account for the detector re¬ 
sponse to jets of hadrons and compared directly 
to predictions from two Monte Carlo (MC) simu¬ 
lations. The simulations considered are the ALP- 
GEN generator 0 , which calculates the lowest 
perturbative order of the strong interaction that 
produces a W + W~ pair of vector bosons and a 
fixed number of partons in the final state (typ¬ 
ically referred to as fixed-order MC), and the 
MC@NLO next-to-leading order (NLO) genera¬ 
tor 0 both interfaced to parton-shower genera¬ 
tors. Fixed-order and NLO simulation are widely 
used methods for modeling electroweak processes 
with multiple jets. The measurement of differ¬ 
ential cross sections as a function of the num¬ 
ber of associated jets is particularly suited to the 
Tevatron because the large rate of top-quark-pair 
production at the LHC makes this analysis much 
more difficult. The production of VF-boson pairs 
was first observed by the CDF experiment using 
Tevatron Run I data [l2j . This process has since 
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been measured by the CDF 0], DO [l3| . AT¬ 
LAS 0, and CMS 0,113 experiments. 

This measurement uses an integrated luminos¬ 
ity of 9.7 fb _1 of pp collision data collected by the 
CDF experiment which comprises the full Teva- 
tron Run II data set. A 12% increase in signal ac¬ 
ceptance is obtained by using events in which one 
or both W bosons decay to a tau lepton that sub¬ 
sequently decays to an electron or muon. Events 
are classified based on the multiplicity and the 
transverse energy of jets, where jets are the prod¬ 
ucts of the parton showering and hadronization of 
high energy partons produced in the initial scat¬ 
tering process. The cross section measurement 
uses an artificial neural network (NN) method to 
distinguish signal and background events. The 
jet multiplicity classifications are zero, one, and 
two or more jets and separate NNs are used for 
each multiplicity class. Events with a single jet 
are further classified according to the energy of 
the jet transverse to the beam line into bins of 
15—25 GeV, 25—45 GeV, and more than 45 GeV. 
The cross sections are extracted via a maximum- 
likelihood fit to the data of a weighted sum of 
the normalized binned NN distributions for sig¬ 
nal and backgrounds, simultaneously over the 
five event classes. 

The CDF experiment consists of a solenoidal 
spectrometer with a silicon tracker and an open¬ 
cell drift chamber surrounded by calorimeters 
and muon detectors 0]. The kinematic prop¬ 
erties of particles and jets are characterized us¬ 
ing the azimuthal angle (f> and the pseudorapid¬ 
ity p = — ln(tan(0/2)), where 9 is the polar angle 
relative to the nominal proton beam axis. Trans¬ 
verse energy, Et , is defined to be Esin 9, where 
E is the energy deposited in pointing-tower ge¬ 
ometry electromagnetic and hadronic calorime¬ 
ters. Transverse momentum, pt-, is the momen¬ 
tum component of a charged particle transverse 
to the beam-line. 

This analysis uses jets, electrons, muons and 
missing transverse energy. Electron and muon 
candidates (hereafter electrons and muons) are 
typically identified using the drift chamber and 
electromagnetic calorimeter or muon chambers, 
respectively. Jet candidates (jets) are measured 
using an iterative cone algorithm that clusters 
signals from calorimeter towers within a cone of 
A R = yj (A rj) 2 + (A <j)) 2 = 0.4. Corrections are 
applied to the calorimeter cluster-energy sums to 
account for the average calorimeter response to 
jets [0. The missing transverse energy vector, 

Ift T , is defined as the opposite of the vector sum 


of the Et of all calorimeter towers, corrected to 
account for the calorimeter response to jets and 
to muons (3j. 

The experimental signature of the signal events 
is two leptons with opposite charge and large If T . 
A number of SM processes result in this observed 
final state and are thus backgrounds. The WZ 
and ZZ processes can yield the same final state 
if one boson decays hadronically, or some decay 
products are not detected. Top quarks decay pre¬ 
dominantly to a b quark and a W boson, which 
can decay leptonically, making ti production a 
dominant background for events with two jets. 
Incorrect measurements of the energy or momen¬ 
tum of leptons or jets result in apparent I)j T . 
This allows the large-production-rate Drell-Yan 
(DY) process, Z /y* —>■ £ + £~, where i refers to 
any charged lepton, to be included in the candi¬ 
date sample. A third source of background arises 
from the misidentification of a final-state parti¬ 
cle, such as in IV+jets and Wy processes, where 
the photon or jet is incorrectly reconstructed as 
a lepton. 

Events are selected from a sample containing 
at least one electron (muon) of Et(pt) > 20 GeV 
and I 77 I < 2.0 that satisfy an online single-lepton 
selection requirement in the trigger [20]. A sec¬ 
ond lepton is required to have opposite charge 
and Et{pt) > 10 GeV. Jets are required to 
have Et > 15 GeV and \r]\ < 2.5. In order 
to reduce the contribution of DY events with 
mismeasured Et(pt), the I$T,rei defined as 

IjjT rei — $t sin ^(^ri ^/j) if the azimuthal 

separation between the I} T an( i the momentum 
vector of the nearest lepton or jet is less than 
Otherwise rei = rfr- rei is required to be 
greater than 25 GeV for same-flavor lepton pairs, 
or 15 GeV for electron-muon events from the 
leptonic decays of r + r— DY events, where the 
DY contribution is small. Requiring the invari¬ 
ant mass of the same-flavor lepton pairs () to 
be outside the Z mass range (between 80 and 99 
GeV/c 2 ) further suppresses Z /y* -» £ + £~ events, 
while requiring one radian of azimuthal separa¬ 
tion between the If T and the dilepton momen¬ 
tum suppresses Z /y* — > t + t~ events. Leptons 
are required to be isolated, by applying the crite¬ 
ria that the sum of Et for calorimeter towers (or 
the sum of the pt of charged particles for cen¬ 
tral muons) in a cone of A R around the lepton 
is less than 10 % of the electron Et (muon pr). 
This helps to purify the lepton sample and reduce 
backgrounds due to misidentified particles, par- 
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ticularly kP±jets and W 7 . The W'y background 
is further reduced by requiring mu greater than 
16 GeV/c 2 . To suppress tt background we re¬ 
ject events with two or more jets in which any 
jet passes a 5-flavored-quark identification algo¬ 
rithm ( 2 lj . 

The geometric and kinematic acceptance for 
the WW signal and the yields of the DY, WZ, ZZ, 
IP 7 and tt backgrounds is determined by simu¬ 
lation. Backgrounds consisting of DY leptons as¬ 
sociated with the production of zero or one jet 
are simulated using the pythia generator |22 |. 
DY background with two or more jets and WW 
samples are simulated with up to three partons 
with the ALPGEN generator. The contributions 
to the WW production from the next-to-next- 
to-leading-order (NNLO) process gg —> W + W~ 
are accounted for by reweighting the WW sim¬ 
ulation as a function of the WW pr distribu¬ 
tion to incorporate the extra contribution pre¬ 
dicted in Ref. [23]. The DY background is stud¬ 
ied in a sample of low Ifi T rel events and a cor¬ 
rection is applied to account for the mismodeling 
of ljj T in the simulation of zero-jet and one-jet 
events Q. The WZ, ZZ, and tt samples are sim¬ 
ulated with the pythia generator [22]. The ex¬ 
pected yields for simulated processes are normal¬ 
ized using cross section calculations performed 
at NNLO for tt [Hj], NLO for WZ, ZZE^ and 
DY [26|, and fixed order for WW [lo|T The 
IP 7 background is simulated using the Baur MC 
generator [27| . but is normalized according to a 
study of same-charge lepton with low mu Pal - 
Samples are generated using the CTEQ5L [28| 
parton distribution functions (PDFs), and in¬ 
terfaced to pythia for parton showering, frag¬ 
mentation, and hadronization. The detector re¬ 
sponse for these processes is modeled with a 
GEANT3-based simulation (29|, and further cor¬ 
rected for trigger, reconstruction, and identifica¬ 
tion efficiencies using samples of W —> ev and 
Z —\ it events [20]. The kF±jets background is 
determined using collision data. The probability 
for a jet to be reconstructed as a lepton is mea¬ 
sured in jet-triggered data samples and applied 
to a kP±jets data sample to estimate the contri¬ 
bution from PP+jets events. The estimated and 
observed yields for signal and background in each 
jet region are shown in Table |U 

We use a NN method to further discrimi¬ 
nate signal from background. Separate NNs 
are trained using the NEUROBAYES [3C| program 
with simulated signal and background events for 


Process Events (best fit) 


Zero jets One jet Two or more jets 


WZ 

19.5 ±3.0 

16.7 ±2.3 

4.26 ±0.81 

ZZ 

13.2 ± 1.9 

4.25 ±0.61 

1.33 ±0.26 

tt 

3.7 ± 1.0 

76 ± 12 

158 ±16 

DY 

150 ± 34 

83 ±21 

20.2 ± 8.6 

Wj 

214 ±27 

44.0 ±6.4 

7.5 ± 1.9 

kP±jets 

685 ± 118 

250 ± 46 

81 ±15 

Total bk 

1086 ± 124 

474 ± 57 

272 ± 26 

IPIP 

963 ± 108 

224 ± 29 

73 ±20 

kPkP+bk 2049 ± 177 

698 ± 73 

345 ± 39 

Data 

2090 

682 

331 

TABLE I: 

Estimated 

and observed 

event yields. 


Event yields and uncertainties are normalized to the 
values returned by fitting signal and background (bk) 
distributions to the data. Uncertainties are due to 
statistical and systematic sources as described in the 
text. 



NN output 


FIG. 1: Estimated and observed combined distribu¬ 
tions of outputs for neutral networks (NNs) trained 
to identify WW events, with a higher NN output 
value indicating a more signal-like event. Event 
yields are normalized to the values returned by fit¬ 
ting signal and background distributions to the data. 


events with zero, one, and two or more jets. In¬ 
puts to the NNs are selected to exploit features 
of the signal and background processes. The in¬ 
puts include the total energy of the event, which 
tends to be modest in WW events compared to 
tt events; the pr of the second highest pt lepton, 
which tends to be low in events where a jet or 
photon is misidentified as a lepton; missing Et 
that is not aligned with a lepton or jet, suggest¬ 
ing that it is the result of neutrino production; 
and, in events with two or more jets, the vec- 
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tor sum of the pr of the leading two jets, which 
tends to be greater if the jets are tt decay prod¬ 
ucts. For events without jets, we use a likelihood 
ratio of the signal probability density over the 
sum of signal and background probability densi¬ 
ties as an input to the NN [1|. The probability 
densities are calculated on an event-by-event ba¬ 
sis for the signal and each background hypoth¬ 
esis using the observed kinematic properties of 
the events in leading-order (LO) matrix-element 
calculations (25j| . The sum of final NN outputs 
for all event classes is shown in Fig. [TJ 

Systematic uncertainties affecting both the 
normalization and shape of the signal and back¬ 
ground NN distributions are assessed for the 
following sources. Uncertainties on acceptance 
originating from lepton-selection and trigger- 
efficiency measurements contribute a 4.3% un¬ 
certainty on all event yields. Acceptance un¬ 
certainties due to potential contributions from 
higher-order effects are evaluated as follows. For 
the WW contribution, the ALPGEN sample is 
reweighted by the pt of the WW system accord¬ 
ing to samples generated with different choices of 
PDF, renormalization, and factorization scales. 
The PDF uncertainty is evaluated with the 
CTEQ61 PDF error set, and ranges from 1.2% 
to 1.8%, increasing as a function of jet multiplic¬ 
ity and jet energy. The scale uncertainty is eval¬ 
uated by doubling or halving the renormaliza¬ 
tion and factorization scales simultaneously and 
ranges from 0.5% to 24%, increasing as a function 
of jet multiplicity and decreasing as a function of 
jet energy. The scale uncertainties are anticor¬ 
related between adjacent bins of jet energy and 
jet multiplicity. Between the zero-jet and one- 
jet bins the anticorrelation is specifically between 
the bin with no jets and one jet but having the 
lowest Et- The effect of the scale uncertainty 
on the shape of the WW NN distributions is also 
evaluated. For tt contribution an uncertainty of 
2.7% is assigned due to QCD effects [|l|]. For 
WZ and ZZ backgrounds, which are simulated at 
LO, a 10% uncertainty is assigned, arising from 
the difference in the observed acceptance of WW 
events generated at LO and NLO with PYTHIA 
and MC@NLO, respectively. Because modeling of 
higher-order amplitudes can affect the extrapo¬ 
lation of the normalization to the predicted W') 
event yield, this uncertainty is also applied to the 
W'y background. A 6.8% uncertainty is also as¬ 
signed to the W'y background due to the photon- 
conversion modeling. An uncertainty of 19% to 
26%, increasing as both a function of jet mul¬ 


tiplicity and jet energy, is assigned to the DY 
background to account for mismodeling of I} T 
based on the differences in acceptance observed 
when varying the I$ T template in the study of 
low rel DY events described above. Uncer¬ 
tainty on jet modeling for simulated backgrounds 
varies from 1.0% to 29%, and are anticorrelated 
between jet multiplicity and jet energy and gen¬ 
erally increases as a function of jet multiplicity. 
For WW and DY processes, the effect on the 
shape of the distribution is also evaluated. For 
the kF+jets background, systematic uncertain¬ 
ties of 20% to 30%, depending on lepton types, 
are determined by calculating the misidentifica- 
tion probabilities at different jet-energy thresh¬ 
olds. Theoretical uncertainties of 6.0% are as¬ 
signed to WZ and ZZ cross sections @ and 4.3% 
to the tt cross section (24|. The measured lumi¬ 
nosity, which is used to normalize the yield of all 
signal and background processes modeled using 
simulations, has an uncertainty of 5.9% [HJ. Ta¬ 
bles including the systematic uncertainties used 
in each jet classification and all correlations and 
anticorrelations are given in Ref. [If]. The sum 
of uncertainties from statistical and systematic 
sources is included for each signal and back¬ 
ground for each jet multiplicity in Tabic HI The 
uncertainties on the background predictions are 
those after the fit described below. 

The differential WW cross section is extracted 
from the NN output shapes incorporating the 
normalizations and systematic uncertainties of 
signal and background in each signal region via 
a binned-maximum-likelihood method. The like¬ 
lihood is formed from the Poisson probabilities 
of observing nt events in the ith bin, in which 
Pi are expected. Systematic uncertainties are in¬ 
cluded as normalization parameters on signal and 
background and subject to constraint by Gaus¬ 
sian terms. The likelihood is given by 


i 


7li — if • 

Mi e 



(1) 


where 


» = Y'**Ifr+r k s e )m<. ( 2 ) 

k c 

For a process k and a systematic effect c, /£ is 
the fractional uncertainty assigned, and S c is a 
Gaussian-constrained floating parameter associ¬ 
ated with c. ( Nk)t is the expected number of 



events from process k in bin i. Systematic un¬ 
certainties that affect the shape of the distribu¬ 
tion are treated as correlated with the appro¬ 
priate rate uncertainties. The parameter au is 
an overall normalization that is fixed to one for 
all processes except WW, for which it is deter¬ 
mined by the fit independently for each anal¬ 
ysis region. The likelihood function is maxi¬ 
mized with respect to the systematic parame¬ 
ters S c and cross section normalizations aww 
simultaneously in all regions. The cross section 
in each region is calculated by multiplying the 
value of aww by the predicted cross sections 
calculated by ALPGEN. The total cross section 
is determined to be cr(pp —> W + W~ + X) = 
14.0 ± 0.6 (stat) 1 ^'q (syst) ± 0.8 (lumi) pb which 
is consistent within one a with the inclusive NLO 
cross section prediction of 11 ± 0.7 pb as cal¬ 
culated by the MC@NLO program and the total 
prediction of the fixed-order program ALPGEN. 
The result is unfolded to the hadronic-jet level 
based on the results of a study of the bin-to- 
bin migration of events due to jet reconstruc¬ 
tion, jet-energy scale, and jet-resolution effects 
as determined in simulated events. The final re¬ 
sult is iteratively corrected to account for differ¬ 
ences in acceptance between the reconstructed 
and true distributions using a Bayesian [33l l34j 
technique. Migrations between jet-multiplicity 
and jet-energy bins are typically of order 10% 
or less. An independent training sample of sim¬ 
ulated events is used to test the unfolding pro¬ 
cess. The correct differential cross sections are 
reproduced stably with a minimal number of it¬ 
erations. The unfolded results are compared to 
ALPGEN, using the CTEQ5L PDFs and inter¬ 
faced to PYTHIA for parton showering with the 
Michelangelo L. Mangano (MLM) matching al¬ 
gorithm 35j, and MC@NLO, usin g th e CTEQ5M 
PDFs and interfaced to herwig j36l. [37j for par- 
ton showering. The measured and predicted dif¬ 
ferential cross sections are shown in Table |TT] and 
Fig. [2] The differential cross section measure¬ 
ments are consistent with both simulation predic¬ 
tions. The largest deviation occurs in the two-jet 
sample being less that two standard deviations. 
NNLO contributions to the qq -A W + W~ cross 
section are not accounted for in the calculations 
and are expected to increase the predicted cross 
section. 


In summary, the WW cross section is measured 
in the dilepton channel both inclusively and dif¬ 
ferentially in jet multiplicity and Et using a 



FIG. 2: Measurement and predictions of a(pp —► 
W + W~ + njets). Values are given inclusively and 
differentially as functions of jet multiplicity and jet- 
transverse energy. Transverse energy ranges are (a) 
15 < Et < 25 GeV, (b) 25 < Et < 45 GeV, and (c) 
Et > 45 GeV. 


neutral-net discriminant and binned-maximum- 
likelihood fit. This is the first measurement of 
the differential cross section for pair production 
of massive-vector bosons. The measured cross 
section, 14.0 ± 0.6 (stat) 1 ^'q (syst) ± 0.8 (lumi) 
pb, and the differential cross sections are con¬ 
sistent with both the NLO and fixed-order pre¬ 
dictions. This result indicates the suitability 
of using either of these theoretical techniques 
to study processes with multiple gauge bosons 
and jets. Processes of this type will be used 
extensively at the LHC to perform searches for 
non-SM physics and to investigate the nature of 
clectroweak-symmetry breaking by studying the 
process of vector-boson scattering. 
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cr(pb) 


Uncertainty (pb) 


o-(pb) 

Jet bin 

Measured Statistical Systematic Luminosity 

ALPGEN MC@NLO 

Inclusive 

14.0 

±0.6 

+ 1.2 
-1.0 

±0.8 

11.3 ± 1.4 11.7 ± 0.9 

Zero jets 

9.57 

±0.40 

+0.82 

-0.68 

±0.56 

8.2 ± 1.0 8.6 ± 0.6 

One jet inclusive 

3.04 

±0.46 

+0.48 

-0.32 

±0.18 

2.43 ± 0.31 2.47 ± 0.18 

One jet, 15 < E T < 25 GeV 

1.47 

±0.17 

+0.13 

-0.09 

±0.09 

1.26 ± 0.16 1.18 ± 0.09 

One jet, 25 < E T < 45 GeV 

1.09 

±0.18 

+0.14 

-0.11 

±0.06 

0.77 ± 0.10 0.79 ± 0.06 

One jet, Ep > 45 GeV 

0.48 

±0.15 

+0.19 

-0.11 

±0.03 

0.40 ±0.05 0.46 ±0.03 

Two or more jets 

1.35 

±0.30 

+0.45 

-0.28 

±0.08 

0.64 ± 0.08 0.61 ± 0.05 


TABLE II: Measurements and predictions of a(pp —► W + W + jets). Values are given inclusively and 
differentially as functions of jet multiplicity and jet-transverse energy. 
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